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The non-axial reflection-asymmetric $32 shape in some transfermium nuclei with N = 150, namely 
246 Om, 748Cf, ?50Fm, and 7°?No are investigated with multidimensional constrained covariant den- 
sity functional theories. By using the density-dependent point coupling covariant density functional 
theory with the parameter set DD-PC1 in the particle-hole channel, it is found that, for the ground 
states of 745Cf and 7°°Fm, the non-axial octupole deformation parameter 632 > 0.03 and the energy 
gain due to the 632 distortion is larger than 300 keV. In 74°Cm and ?°?No, shallow $32 minima are 
found. The occurrence of the non-axial octupole {32 correlations is mainly from a pair of neutron 
orbitals [734]9/2 (vji5/2) and [622]5/2 (vgo/2) which are close to the neutron Fermi surface and a 
pair of proton orbitals [521]3/2 (af7/2) and [633]7/2 (i13/2) which are close to the proton Fermi 
surface. The dependence of the non-axial octupole effects on the form of energy density functional 


and on the parameter set is also studied. 
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The majority of observed nuclear shapes is of 
spheroidal form which is usually described by 629. The 
existence of the nonaxial-quadrupole (triaxial) deforma- 
tion B22 [1-3] and the axial octupole deformation 830 [4] 
in atomic nuclei have also been confirmed both experi- 
mentally and theoretically. However, there is no a pri- 
ori reason to neglect the nonaxial-octupole deformations, 
especially the 632 deformation [5-7]. The pure (32 de- 
formation has a tetrahedral symmetry with a symme- 
try group TP. The existence of non-trivial irreducible 
representations of this group makes it possible for a nu- 
clei to have large energy gaps in their single-particle lev- 
els, thus increasing it’s stability [8]. It has been antici- 
pated that 832 deformation occurs in the ground states 
of some nuclei with special combinations of the neu- 
tron and proton numbers [7-9]. Recently, lots of theo- 
retical studies focus on this nuclear shape, either from 
the TP-symmetric single particle spectra [7, 9-11] or 
from various nuclear models including the macroscopic- 
microscopic model [9, 11-13], the Skyrme Hartree-Fock 
(SHF), SHF+BCS, or Skyrme Hartree-Fock-Bogoliubov 
models [11-18], and the Reflection Asymmetric Shell 
Model (RASM) [19, 20]. In particular, Dudek et al. 
predicted that a negative-parity band in !°°Gd is a fa- 
vorable candidate to manifest tetrahedral symmetry [13] 
which has stimulated several interesting experimental 
studies [21, 22]. 

Nowadays the study of nuclei with Z ~ 100 becomes 
more and more important because it not only reveals the 
structure of these nuclei themselves but also provides sig- 
nificant information for superheavy nuclei [23-26]. One 
of the relevant and interesting topics is how to explain the 
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low-lying 27 states in some N = 150 even-even nuclei. 
In these nuclei, the bandhead energy E(27) of the lowest 
27 bands is very low [27]. It is well accepted that the 
octupole correlation is responsible for it. For example, 
a quasiparticle phonon model with octupole correlations 
included was used to explain the excitation energy of the 
27 state of the isotones with N = 150 and the octupole 
correlation is mainly due to the neutron two-quasiparticle 
configuration 9/27 [724] @ 5/2*[622] and proton two- 
quasiparticle configurations 9/2*[633] @ 5/27[521] or 
7/2*[633] 9 3/27[521] [28]. In all these configurations, 
the third components of the angular momenta of like- 
quasiparticles K differ by 2, i.e., Y32-correlation should 
play an important role. In Ref. [20], Chen et al. pro- 
posed that the non-axial octupole Y32-correlation results 
in the experimentally observed low-energy 2~ bands in 
the N = 150 isotones and the RASM calculations re- 
produces well the experimental observables of these 27 
bands. It was also predicted that the strong nonaxial- 
octupole effect may persist up to the element 108 [20] 
and play a crucial role in determining the shell stabil- 
ity in even heavier nuclei [29]. Pronounced shell effects 
were also found for N = 16, 40, and 110 when combining 
nonaxial octupole deformations [30]. 

In this Brief Report we present a microscopic and self- 
consistent study of the Y32 effects in the N = 150 isotones 
under the framework of the covariant density functional 
theory (CDFT) [31-36]. We use the recently developed 
multi-dimensional constraint (MDC) CDFT [87, 38] in 
which not only the axial [39, 40] but also the reflection 
symmetries [41] are broken. For the parametrization of 
the nuclear shape, the conventional ansatz in mean-field 
calculations 


At 


Pru = SAR 


(Qau); (1) 


201708.00248v1 


chinaXiv 


TABLE I. The quadrupole deformation G20, non-axial oc- 
tupole deformation 632, and hexadecapole deformation (40 
together with binding energies E.. for the ground states of 
N = 150 nuclei calculated with DD-PC1. Eaepth denotes the 
energy difference between the ground state and the point with 
832 = 0 in the potential energy curves shown in Fig. 1. Eexp 
denotes experimental binding energies taken from [45]. All 
energies are in MeV. 


Nucleus 620 (32 = Bao Eval Edepth Eexp 

z6 Cm 0.296 0.020 0.126 —1847.932 0.034 —1847.819 
248Cf 0.299 0.037 0.111 —1857.853 0.351 —1857.776 
250Fm 0.293 0.034 0.097 —1865.843 0.328 —1865.520 
252No 0.293 0.025 0.083 —1872.133 0.104 —1871.305 


is adopted, where Q),, = TV is the mass multipole op- 
erator. The nuclear shape is assumed to be invariant un- 
der the reversion of x and y axes in MDC-CDFT, i.e., the 
intrinsic symmetry group is V4 and all shape degrees of 
freedom fà, deformations with even u are allowed. The 
CDFT functional can be one of the following four forms: 
the meson exchange or point-coupling nucleon interac- 
tions combined with the non-linear or density-dependent 
couplings [37, 38]. In the present work, the pairing effect 
is taken into account by using the BCS approximation 
with a separable pairing force [42—44]. 


First we study the ground states of nuclei 746Cm, 
2480 f, 250Fm, and 2°?No by using the density-dependent 
point coupling covariant density functional theory with 
the parameter set DD-PC1 in the particle-hole chan- 
nel [46]. The quadrupole deformations 820, non-axial oc- 
tupole deformation (32, and hexadecapole deformations 
Bao together with binding energies Eea; are listed in Ta- 
ble I. The non-axial quadrupole and axial octupole defor- 
mation parameters are all zero for these nuclei. As is seen 
in Table I, the calculated binding energies agree very well 
with the data. All of these four nuclei are well deformed 
with the quadrupole deformation parameter 20 ~ 0.3 
and the hexadecapole deformations G49 œ% 0.1. Super- 
posed with large quadrupole and hexadecapole deforma- 
tions, finite values are obtained for the non-axial octupole 
deformation (32 for these N = 150 isotones. 


In order to see more clearly the development of non- 
axial octupole (32 deformation along the N = 150 iso- 
tonic chain, we preform one-dimensional constrained cal- 
culations and obtain potential energy curves, i.e., the to- 
tal binding energy as a function of (632. At each point of 
a potential energy curve, the energy is minimized auto- 
matically with respect to other shape degrees of freedom 
such as 829, 822, 830, and 640, etc. In Fig. 1, we show 
potential energy curves for these N = 150 isotones. For 
246Cm, the ground state deformation 632 = 0.020 (see 
Fig. 1 as well as Table I). The potential energy curve 
is rather flat around the minimum. We denote the en- 
ergy difference between the ground state and the point 
with 832 = 0 by Eaepth which measures the energy gain 
with respect to the 632 distortion. For 74°Cm, Egeptn is 


E (MeV) 
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FIG. 1. (Color online) The binding energy E (relative to the 
ground state) for N = 150 isotones ?“°Cm (dashed line), °S Cf 
(solid line), 7°°Fm (long-dashed line), and 7°*No (dotted line) 
as a function of the non-axial octupole deformation parameter 


P32. 


only 34 keV. For 748Cf, 2°°Fm and 2°?No, the minima lo- 
cate at 832 = 0.037, 0.034, and 0.025, respectively. The 
corresponding energy gain Eaepth = 0.351, 0.328, and 
0.104 MeV. It may be hard to conclude that these nuclei 
have static non-axial octupole deformations from these 
results because the potential energy curve is flat around 
the minimum and Egepth is small. However, the present 
calculations at least indicate a strong Y32-correlation in 
these nuclei. Both the non-axial octupole parameter 832 
and the energy gain Egepth reach maximal values at 748Cf 
in the four nuclei along the N = 150 isotonic chain. This 
is consistent with the analysis given in Refs. [20, 28] and 
the experimental observation that in ?4°Cf, the 27 state 
is the lowest among these nuclei. 

The triaxial octupole Y32 effects stem from the cou- 
pling between pairs of single-particle orbits with Aj = 
Al = 3 and AK = 2 where j and / are the total and 
orbit angular momentum of single particles respectively 
and K the projection of j on the z axis. If the Fermi 
surface of a nucleus lies close to a pair of such orbitals 
and these two orbitals are nearly degenerate, a strong 
non-axial octupole 832 effect is expected. 

In Fig. 2, we show the proton and neutron single- 
particle levels near the Fermi surface for ?4°Cf as func- 
tions of quadrupole deformation $29 on the left side and 
of 632 with 329 fixed at 0.3 on the right side. In Fig. 2(a), 
one finds a strong spherical shell closure at Z = 92. As 
discussed in Ref. [47], this shell closure is a spurious one 
and it is commonly predicted in relativistic mean field 
calculations. The spherical proton orbitals 72 f7/2 and 
m1i13/2 are very close to each other and this near degen- 
eracy results in octupole correlations. The two proton 
levels, [521]3/2 originating from 2 f7/2 and [633]7/2 orig- 
inating from 1713/2, satisfying the Aj = Al = 3 and 
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FIG. 2. (Color online) The single-particle levels near the 


Fermi surface for (a) protons and (b) neutrons of *48Cf as 
functions of quadrupole deformation {29 (left side) and of B32 
with b20 fixed at 0.3 (right side). 


AK = 2 condition, are very close to each other when 820 
varies from 0 to 0.3. Therefore the non-axial octupole Y32 
develops and with (32 increasing from zero, an energy 
gap appears at Z = 98. The appearance of the spurious 
shell closure at Z = 92 is mainly due to the lowering 
of w1hg/z orbital which also results in that the proton 
[514]7/2 level lies in between [521]3/2 and [633]7/2, thus 
making the gap at Z = 98 smaller and weakening the 
Yz2 correlation. Similarly, the spherical neutron orbitals 
V2g9/2 and v1ji5/2 are very close to each other. The 
two neutron levels, [734]9/2 originating from 1152 and 
[622]5/2 originating from 2g9/2, are also close to each 
other and they just lie above and below the Fermi sur- 
face. This leads to the development of a gap at N = 150 
with $32 increasing. Therefore it is clear that the Y32 
correlation in N = 150 isotones is from both protons and 
neutrons and for 748Cf the correlation is the most pro- 
nounced. The gap around N = 150 is larger than that 
around Z = 98, which may indicate that the non-axial 


TABLE II. The quadrupole deformation (20, octupole de- 
formation (32, hexadecapole deformation 649 together with 
binding energies Æ for the ground states of N = 150 nuclei 
calculated with parameter sets PC-PK1, DD-ME1, DD-ME2 
and DD-PC1. Eaepth denotes the energy difference between 
the ground states and the point constrained (32 to zero. All 
energies are in MeV. 


Nucleus Parameters $20 (32 (ao E Edepth 
26Cm  PC-PK1 0.304 0.0 0.129 —1846.273 — 
DD-ME1 0.297 0.0 0.126 —1847.205 — 
DD-ME2 0.296 0.0 0.125 —1846.662 — 
DD-PC1 0.296 0.020 0.126 —1847.932 0.034 
2748 Cf PC-PK1 0.306 0.031 0.113 —1857.327 0.139 
DD-ME1 0.299 0.035 0.107 —1857.718 0.278 
DD-ME2 0.300 0.039 0.105 —1857.088 0.474 
DD-PC1 0.299 0.037 0.111 —1857.853 0.351 
2>0Rm  PC-PK1 0.302 0.029 0.098 —1865.301 0.128 
DD-ME1 0.293 0.033 0.093 —1866.425 0.305 
DD-ME2 0.292 0.034 0.091 —1865.505 0.416 
DD-PC1 0.293 0.034 0.097 —1865.843 0.328 
226 PC-PK1 0.300 0.013 0.085 —1870.414 0.001 
DD-ME1 0.293 0.023 0.080 —1873.412 0.089 
DD-ME2 0.292 0.024 0.078 —1872.235 0.125 
DD-PC1 0.293 0.025 0.083 —1872.133 0.104 


octupole effect originating from neutrons is larger than 
that from protons. If there is not a spurious shell gap 
at Z = 92, one can expect more pronounced non-axial 
octupole correlations from protons. 


In order to examine the dependence of our results on 
the functional form and on the effective interaction, we 
also studied 74°Cm, 748Cf, 2°°Fm, and ?52No with other 
covariant density functionals and several other typical 
parameter sets, including PC-PK1 [48], DD-ME1 [49], 
and DD-ME2 [50]. The results are listed in Table IT. 


Roughly speaking, the results are similar with differ- 
ent parameter sets: $22 and {39 vanish in all cases and 
for a specific nucleus, the calculated 629 and {49 and the 
binding energy E with different parametrizations agree 
with each other. For 74°Cm, three parameter sets PC- 
PK1, DD-ME1, and DD-ME2 predict zero $32. With 
DD-PC1, as discussed before, a rather shallow poten- 
tial energy curve with a minimum at (32 = 0.020 is 
obtained. This indicates that the octupole correlation 
in 746Cm is very weak. Both from the value of 632 
and from the energy gain Egepth due to the (32 distor- 
tion, it is found that the density-dependent functionals 
gives stronger Y32 effects than the functional with non- 
linear self coupling does. Among the density-dependent 
functionals, the meson exchange nucleon interaction with 
the parameter set DD-ME2 gives the largest energy gain 
Eaepth for 48Cf£, 25°Fm, and 752No. The evolution of the 
non-axial octupole (32 effects along the N = 150 isotonic 
chain is almost independent of the form of energy den- 
sity functional and the parameter set: The Y32 effects is 
the strongest in ?4°Cf except for that DD-ME1 gives a 
smaller energy gain for 748Cf than for 7°°Fm. 
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In summary, we studied the non-axial octupole Y32 
effects in N = 150 isotones 74°Cm, 748Cf, 75°Fm, and 
252No by using multi-dimensional constraint covariant 
density functional theories. Due to the interaction be- 
tween a pair of neutron orbitals, [734]9/2 originating from 
Vji5/2 and [622]5/2 originating from vgg/2, and that of 
a pair of proton orbitals, [521]3/2 originating from r f7/2 
and [633]7/2 originating from 73/2, rather strong non- 
axial octupole Y32 effects have been found for 748Cf and 
250Fm which are both well deformed with large axial- 
quadrupole deformations, 329 ~ 0.3. For 74°Cm and 
252No, a shallow minima develops along the 832 deforma- 
tion degree of freedom. The evolution of the non-axial 
octupole (32 effect along the N = 150 isotonic chain is 
not very sensitive to the form of the energy density func- 
tional and the parameter set we used. Finally we note 
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